Storage and Retrieval of a Microwave Field in a Spin Ensemble 
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We report the storage and retrieval of a small microwave field from a superconducting resonator 
into collective excitations of a spin ensemble. The spins are nitrogen-vacancy centers in a diamond 
crystal. The storage time of the order of 30 ns is hmited by inhomogeneous broadening of the spin 
ensemble. 



Superconducting qubits are promising candidates for 
quantum information processing; however their coher- 
ence times [1_ cannot yet compete with those of micro- 
scopic systems such as atoms |2] , electron or nuclear spins 

. Hybrid quantum circuit architectures have thus been 
proposed [4-8_, in which microscopic systems would be 
used as quantum memory for superconducting qubits. 
Whereas the coupling of one individual atom or spin to 
a superconducting circuit is usually too weak, the cou- 
pling constant of an ensemble of N such systems is en- 
hanced by \/N, allowing to reach the strong coupling 
regime requested for quantum information applications. 
Proposals for spin-ensemble based hybrid quantum cir- 
cuits often consist of a superconducting resonator used 
as a quantum bus between the ensemble and the su- 
perconducting qubit. On the experimental side [9-12J, 
strong coupling between a spin ensemble and a supercon- 
ducting resonator has up to now been demonstrated only 
spectroscopically. Here we report the first time-domain 
measurements of the coherent storage and retrieval of a 
classical microwave field jT^ of about 100 photons from a 
superconducting resonator into collective excitations of a 
spin ensemble consisting of negatively-charged nitrogen- 
vacancy centers in diamond (NV centers) , an important 
step towards a spin-based hybrid quantum circuit archi- 
tecture. 

The experiment relies on the fact that the interaction 
between the electromagnetic field in the cavity mode and 
the spin ensemble involves only one collective spin vari- 
able, which behaves as a harmonic oscillator in the limit 
of low excitation energy [8|. This effective spin oscil- 
lator is magnetically coupled with a collective coupling 
constant (/ens to the superconducting resonator whose fre- 
quency can be tuned. When two such coupled harmonic 
oscillators are suddenly put into resonance, they coher- 
ently exchange energy with a period n/gcns- We observe 
this dynamics by measuring the amplitude of the mi- 
crowave field leaking out of the resonator after its in- 
teraction with the spins, which is found to oscillate as 
a function of the interaction time. This storage-retrieval 
cycle is however damped in a relatively short time, which 
as we show is limited by the inhomogeneous broadening 
of the NV center ensemble. Quantitative agreement with 
recent theoretical work [141 [T5] is obtained for a consis- 



tent set of data, covering spectroscopic as well as time- 
domain measurements. 

In our experiment (see Figllk and Supplementary 
Material), an ensemble of ^ 10 NV centers in a di- 
amond crystal are magnetically coupled to a supercon- 
ducting resonator. A static magnetic field Snv — 1-7 mT 
is applied to the spins, parallel to the chip surface, along 
the [1, 0, 0] crystal axis within a few degrees (see Figjljo). 
With this orientation, the four possible NV center crys- 
talline orientations all make approximately the same an- 
gle 9 ~ 55° with -Bnv so that their resonance frequencies 
(see Fig|lj3) are approximately equal. The resonator in- 
cludes a four-SQUID array so that its frequency uj^{^) 
can be tuned by changing the fiux $ through the SQUIDs 
loop [16 . Such a frequency tuning can be done on 
a nanosecond timescale |17j. For time-domain exper- 
iments, we measure the amplitude A{t) of microwave 
pulses transmitted through the resonator using homo- 
dyne detection followed by sampling and averaging. 

We first characterize the resonator-spins system by 
measuring the resonator transmission while scanning $. 
As shown in Figjlji, two vacuum Rabi splittings are ob- 
served when the resonator frequency matches either of 
the two NV center transitions w_(resp. a;+) from ground 
state TOg = to TOs = —1 (resp. tos = +1). A fit of 
these data using a coupled oscillators model [S] yields 
(;ons/27r = 10.6 MHz for the lower frequency anticrossing 
occurring at uj-/2tt = 2.85 GHz, on which we will focus 
in the following. The transmission close to the middle 
of the anticrossing is shown in Fig. [T^ where two well- 
resolved polaritonic peaks can be seen, an indication that 
the two oscillators are in the strong coupling regime and 
that the coherent exchange of excitations between the 
resonator and the spin ensemble can be observed in the 
time-domain. 

To demonstrate such dynamics, the experiment pro- 
ceeds as follows (see Fig|2|i,). At the beginning of each 
experimental sequence (averaged out typically 5 • 10^ 
times) , the resonator is biased at a fiux $off such that 
its frequency a;r(<i>oFF)/27r — 2.88 GHz is out of reso- 
nance with the spins. A microwave pulse at frequency 
<^r(*i'oFF) of duration 2/is (much longer than the res- 
onator damping time Tcav = Q/i^r ^ 100 ns) estab- 
lishes a steady-state coherent field of amplitude a inside 
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Figure 1: (a) Setup description. The NV centers ensemble is 
magnetically coupled to the coplanar resonator containing a 
SQUID. The flux $ through the SQUID loop can be tuned on 
a nanosecond scale by applying current pulses to an on-chip 
antenna. A magnetic field Bnv = 1.7 mT is applied to the 
spins, parallel to the sample and to the [1, 0, 0] crystal axis (b). 
The amplitude A{t) of microwave pulses transmitted through 
the resonator is measured by homodyne detection at room- 
temperature after amplification, (c) Simplified energy level 
scheme of a NV center, (d) Resonator transmission \S2i \ (ijj) 
at -Bnv = 1.7 mT as a function of $ (in units of the supercon- 
ducting flux quantum <I>o = h/2e) showing two vacuum Rabi 
splittings. The solid line is a flt to the data using the coupled 
oscillators model 0. (e) Vacuum Rabi splitting close to res- 
onance with UJ-. Black line is experimental data, red line is 
theory as explained in the text (rescaled in amplitude to flt 
the data), assuming a spins-resonator detuning of 0.5 MHz. 

the resonator (with \a\^ ^ 120 photons). Right after 
the microwave pulse is switched off (t ~ 0), the res- 
onator frequency is brought close to a;_ by a flux pulse 
of amplitude A$ (see Supplementary Material) and du- 
ration T, during which the resonator and the spin en- 
semble may exchange energy causing the intraresonator 
field amplitude \a{t)\ to oscillate. After the flux pulse, 
the only evolution of the field is an exponential decay 
\a{t > r)| = |a(T)|exp(-(t - t)/2T^.^^). Measuring the 
amplitude of the exponentially damped microwave signal 
that leaks out of the cavity therefore directly yields |a(r)| 
and reveals the coupled resonator-spins dynamics. 

Figure [2)3 sh ows the results obtained for a flux pulse 
amplitude A$ that puts the resonator in resonance with 
the spins at u-. The two-dimensional plot shows the 
measured microwave output amplitude A{t) for increas- 
ing T. The curve A{t) for the shortest flux pulse t = 5 ns 
(see Fig. [2j:) showing the microwave field decay after the 
pulse is switched off is only approximately exponential, 
due to a slight nonlinearity of the resonator caused by the 
presence of the SQUID [18J, which we neglect in the fol- 
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Figure 2: (a) Sketch of an experimental sequence (b) Output 
amplitude A(t) measured for increasing fiux pulse duration r, 
for a fiux pulse amplitude such that a;, ($oFF + A$) = uj-. 
The dashed orange line indicates t = r. The two white 
dashed lines indicate the time window on which averaging 
is performed to compute a(r) (see text), (c) Amplitude 
A{t) for gansT ~ 0, TT, 27r, Stt. (d) Normalized amplitude 
a(r) — A{t + toff) (see text). Red dots are experimental data, 
continuous blue line is theory as explained in the text, (in- 
set) Sketch of the physical process leading to damping of the 
oscillations because of coupling of the superradiant state to 
the bath of subradiant states, (e) ODMR signal measured at 
room-temperature (red dots), compared to the spin distribu- 
tion p{ijj) (see text) used in the theoretical calculations (blue 
line). 



lowing. For increasing r, the curves A{t) clearly display 
two parts : from t = to t = t rapid transient oscilla- 
tions are observed, which are unfortunately difficult to 
quantitatively interpret due to the flnite bandwidth of 
our homodyne detection setup. After t — t, A{t) shows 
a decay similar to the one observed for the shortest pulse 
but with an amplitude that oscillates with r with a pe- 
riod T ~ 50 ns, in quantitative agreement with the cou- 
pling strength estimated from the vacuum Rabi splitting 
T ~ 7r/(7cns- This establishes that the observed oscillation 
is indeed due to cycles where the microwave fleld is flrst 
stored into a collective spin excitation (0 < gens''" < ^r) 
then retrieved (tt < QensT < Stt). Curves corresponding 
to various steps of the cycle are also shown in Fig. [2]:. 
Note that by using a homodyne detection scheme we can 
only measure a fleld that has a well-deflned phase relation 
with the local oscillator, implying that phase coherence 
is indeed preserved even after several (storage, retrieval) 
cycles. A full quantum state tomography would however 
be needed in order to assess the fldelity of the fleld stor- 
age at the quantum level, which is beyond the scope of 
the present work. 

In order to obtain more quantitative insight into these 
data, we compute for each r the quantity a(T) = 
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A{t + toff) where ioff = 140 ns is an offset that was cho- 
sen to avoid taking into account transients and where an 
additional time-averaging is performed on a 20 ns window 
(see Fig. |2 a). Since a(T) is directly proportional to |a(r)| 
it is a good quantity to investigate the spin-resonator dy- 
namics. It is plotted in Fig|2}i for the previous flux pulse 
parameters. Several cycles are observed before the oscil- 
lation amplitude and offset decay simultaneously to zero. 
The fact that the field amplitude decays to zero within a 
hundred of nanoseconds might come as a surprise given 
the fact that the energy damping time of individual spins 
is in the millisecond range at room-temperature and is 
expected to be well above one second at low tempera- 
ture [111. The only energy dissipation should thus be 
caused by damping of the resonator, which occurs on a 
significantly longer timescale. 

Understanding this decay requires to take into account 
the inhomogeneous broadening of the spin ensemble. Fol- 
lowing [Xr, T5^, we model each spin in the Holstein- 
Primakoff approximation (valid in our experiment which 
is deep in the low-excitation regime) by a harmonic oscil- 
lator of frequency wj with annihilation (creation) opera- 
tor bj (&|), the resonator field being described by the anni- 
hilation (creation) operator a (a^). The system Hamilto- 

nian is then H/h = u,-{^)a^ a+J^'^jblbj+J^ ffj (^lo+^jO^^)! 
gj being the coupling constant of spin j with the res- 
onator. Turning to a new orthogonal basis consisting of 
the superradiant spin-wave mode b = (l/gons) 5j^j °f 
frequency a;_ and TV — 1 orthogonal dark modes dj of 

frequency (Dj, with gons — (X) ISjl^)^^'^' the Hamiltonian 
becomes H/h = CLir($)a^a + uj^b^b + gcns{o.^b + h.c) + 

X^Wjdj^j + E7j(^^c?j + h.c) with 7j - A/y/N. Here A 
is the variance of the distribution, i.e. the inhomo- 



geneous spins linewidth. Note that in this expression we 
have neglected terms coupling the dark modes together. 
From the previous expression we see that the superradi- 
ant spin-wave mode is the only one coupled to the res- 
onator field (with constant gons), but that in the presence 
of inhomogeneous broadening (A > 0) it is also coupled 
to the — 1 dark modes. In the large- limit, these 
dark modes act as a bath into which energy is damped 
at a rate 7inh, explaining the fast decay in Fig|2jl. In- 
troducing the normalized spin density distribution p{uj) 
(centered around instead of a;_), one can show that 
7i„h ^ 27rp(5ons)5ens CS] in the ^cns > A limit. An 
important consequence of this formula is that increas- 
ing ^cns (by e.g. increasing the spins density) actually 
helps reducing 7ijjh provided p{uj) decays faster than oj"^ 
|14l [T5] . This implies that a strong coupling to the cav- 
ity protects the system against the consequences of in- 
homogeneous broadening, an effect called cavity protec- 
tion in [15 . In our experiment however we are not in 
the gens 3> A limit (see below) and these approximate 
formula are not valid. We thus resort to explicit an- 
alytical formula obtained using input-output relations 
on H [TS], which allow us to calculate both the res- 
onator transmission t(uj) = K/[2i(uj — uj,-) — n — 2iW{uj)] 
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Figure 3: Coherent field exchange between the resonator and 
the spins for varying flux pulse ampHtude A$. Normahzed 
ampfitude a(r) (see text) is plotted versus r and A$. Left 
graph is experimental data, right graph is theory. The white 
dashed line indicates A<E>. 



with W{uj) — g'^^^^ p{uj')duj' /[uj — uj'] and the time- 
domain signal as the modulus of the inverse Fourier- 
Laplace transform of t(uj) (see |15| and Supplementary 
Material) . 

To estimate these quantities in our experiment, we 
need to know the spin density distribution p{uj). In 
NV centers ensembles, the inhomogeneous linewidth is 
caused by dipolar interactions with neighboring spins, ei- 
ther paramagnetic nitrogen impurities (so-called PI cen- 
ters) that were not converted into NV centers during 
the sample processing, or ^^C nuclei that are present 
in natural abundance (1.1%) [T^]. Note that neigh- 
boring NV centers do not contribute to the inhomoge- 
neous linewidth because at the temperature of our ex- 
periment they are frozen in the mg — state. In our 
sample, a FWHM hnewidth A/27r 7 MHz was mea- 
sured by optically-detected magnetic resonance (ODMR) 
at room-temperature (see Figj2^), compatible with the 
linewidth expected from its nominal PI centers concen- 
tration ~ 100 ppm |19J. As shown in |14l [T5]. quantita- 
tive predictions for the system dynamics require to know 
not only the overall linewidth but also the detailed shape 
of p{uj). In particular Gaussian and Lorentzian distribu- 
tions yield very different results. Spin ensembles inhomo- 
geneously broadened by dipolar interactions are expected 
to show a Lorentzian lineshape [20 with a cutoff; however 
complications in our experiment arise due to hyperfine 
coupling with the ^'^N atom nuclear spin, and to a possi- 
ble misalignment of i?NV with the [1, 0, 0] crystalline axis 
causing the four NV orientations to have slightly differ- 
ent frequencies. As a result, we assume a phenomenolog- 
ical lineshape for our spin ensemble adjusted for reaching 
good agreement both with spectral and time-domain ex- 
perimental data, as a convolution of a Gaussian and a 
Lorentzian profile with respective widths a and 7. The 
parameters chosen in the following {a/2TT — 5.12 MHz 
and 7/27r = 1 MHz) yield a lineshape p{uj) compatible 
with ODMR data although slightly broader, as can be 
seen in Figj2^ (see Supplementary Material) . Using such 
distribution and the formulas above, we obtain quantita- 
tive agreement for spectroscopic (see Fig. ^) as well as 
time-domain (see Fig. kH) measurements. 
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Figure 4: Ramsey fringes experiment. The normalized mi- 
crowave amplitude a{At) after the pulse sequence shown in 
inset is plotted. Red dots are experimental data, black con- 
tinuous line is theory. 

We also study the dependence of the microwave field 
exchange on the resonator-spins detuning by measuring 
a(T) for various A<i> (see Figjsj). Out of resonance, the 
oscillations amplitude is reduced and their frequency in- 
creases as expected for two coupled oscillators. The 
asymmetry in the data between the A$ < A$ and 
A$ > A$ sides is an artefact of the nonlinear depen- 
dence of Wr on A$ and of a residual hybridization of the 
resonator with the spins caused by the finite initial de- 
tuning [uJr{^oFF) — / 9 ^ 3, and is well reproduced 
by theory. The largest discrepancy is observed for pulse 
amplitudes A4> ^ 90m<i>o, where additional features are 
clearly seen in the experiment. We attribute them to a 
small density of NV centers having a ^^C among their 
closest neighbor, known to shift the electron spin fre- 
quency by ±65 MHz due to the hyperfine interaction, 
and that was not taken into account in the calculation. 

We finally perform a Ramsey-like experiment in order 
to quantify the time during which the microwave field can 
be stored in the spin ensemble. Starting as before with 
a steady-state microwave field in the resonator biased at 
^ofFj the resonator is tuned in resonance with the spins 
at < = for a 7r/2 pulse {gcnsT — '"'/S), after which it is de- 
tuned by Auj/2tt = 30 MHz during At, then tuned back 
to resonance for a second 7r/2 pulse, and finally tuned 
back at Wr($oFF), after which a{At) = A{2t + At + tos) 
is measured. Due to the beating between the effective 



spin oscillator and the resonator, oscillations are observed 
in a{At) at frequency Aw as seen in Fig ji] The oscilla- 
tions amplitude decay time ~ 30 ns gives the storage 
time, on the order of ^ 2/ A as expected. The full cal- 
culation, performed as explained above, yields again rea- 
sonable agreement with the measurements. 

As shown by the previous results, inhomogeneous 
broadening appears as a serious obstacle to the successful 
implementation of hybrid quantum circuits based on spin 
ensembles. One first obvious solution is to obtain samples 
with narrower inhomogeneous linewidths. In the specific 
case of NV centers, this requires a very efficient conver- 
sion rate 77 of PI centers initially present in the diamond 
crystal in NV centers. Improvement by one order of mag- 
nitude compared to the sample used in this work (where 

is a few percent) is within reach. This would reduce the 
inhomogeneous linewidth by the same amount, at which 
point inhomogeneous broadening would be dominated by 
the NV centers hyperfine structure (A/27r ~ 4 MHz for 
^'^'N). Another solution is to increase the NV concentra- 
tion even further to reach the regime gens ^ A and rely 
on the cavity protection effect. Finally, it should be pos- 
sible, by applying refocusing techniques used for atomic 
ensemble-based quantum memories [2T], to recover the 
quantum information lost in the dark states into the su- 
perradiant state and then into the resonator. This would 
be the best way to fully benefit from the spin superior 
coherence properties. 

In conclusion we have demonstrated the coherent stor- 
age and retrieval of a classical microwave field of 100 
photons from a superconducting resonator into collec- 
tive excitations of an ensemble of NV centers and dis- 
cussed the implications of inhomogeneous broadening. 
This work opens the way to the operation of a quantum 
memory in which a single microwave photon would be 
coherently stored and retrieved from a superconducting 
qubit into the spin ensemble. 
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Supplementary Material 

I. MEASUREMENT SETUP 

Figure |5ja) shows the full configuration of the room 
temperature apparatus, and the wiring inside the 
fridge. The fast current pulse is generated by an arbi- 
trary function generator AFG3251 (Tektronix) and input 
into the fridge from the port "Flux in" through a 10 dB 
attenuator. After being filtered by low pass filters (LPFl 
and LPF2) at each temperature stage and further attenu- 
ated by 20 dB at 4 K stage, the flux pulse passes through 
the flux line [see the solid arrows in Fig. |6ja) and (b)], 
and goes back to room temperature through low-pass fll- 
ters (LPF3 and LPF4) and a 20 dB attenuator. The end 
of the fast flux circuit ("Flux out") is terminated by a 
broadband 50 il. 

The microwave pulse is generated by mixing contin- 
uous microwave and a fast DC voltage pulse using an 
IQ-mixer[T]. The microwave pulse is input to the port 
"Reso in" after being attenuated by a tunable attenua- 
tor, and reaches to the input port of the chip through 
attenuators at each temperature stage and a band-pass 
pass fllter (BPFl) at 40 mK. After being flltered by a 
band pass fllter (BPF2) and isolated by two circulators, 
the transmitted signal is amplifled by a cryogenic HEMT 
(high electron mobility transistor) amplifler having noise 
temperature ~ 4 K. The amplifled signal is then further 
amplifled at room temperature by two microwave am- 
pliflers, and demodulated into in-(I) and quadrature-(Q) 
phase components using an IQ-demodulator. The de- 
modulated signal is sampled at 500 MS/s by a fast data 
acquisition card Acqiris CC103 (Agilent) after being fur- 
ther amplifled by two ampliflers. 

The resonator chip flxed onto a printed circuit board is 
mounted inside a sample box made of copper [Fig. [sj^b)]. 
A DC magnetic fleld B^v is applied parallel to the chip 
by passing a DC current through an outer superconduct- 
ing coil. The sample box and the coil are surrounded 
by two magnetic shieldings consisting of permalloy tape 
VC6025X (VacuumSchmelze) and a lead cylinder [not 
shown in Fig. [Sjb)]. 



II. FLUX PULSE 

Here we describe a few technical details on the flux 
pulse used to tune the resonator frequency at the 
nanosecond scale. As shown in Fig. [6] the pulse is 
applied using an arbitrary function generator AFG3252 
(Tektronix) that sends a current pulse through an on-chip 
line. The risetime of the voltage pulse is 2 ns. Passing 
current through the on-chip antenna generates screen- 
ing currents through the ground planes, as shown by the 
dashed arrows in Fig. [6]jb), of our superconducting cir- 
cuit with very low damping times (of the order of the 
ms) that in turn affect the flux through the SQUID. As 



a result, the SQUID offset flux bias point was found to 
depend on the time integral of the flux pulse. In addi- 
tion, the flux applied to the SQUID loop by sending a 
fast current pulse through the flux line was strongly re- 
duced (by a factor ^ 50) compared to the flux applied 
with a dc current of the same amount. We measure this 
current-to-flux transfer function with calibration exper- 
iments, allowing us to convert a voltage pulse into flux 
as in flgure 3 of the main text. In order for the bias 
point to be independent of the amplitude and duration 
of the pulse, we had to add, at the end of each sequence, 
a compensation flux pulse opposite to the flrst one [see 
Fig. [6|^c)]. Such a compensation pulse has strictly no 
incidence on the experiment outcome since it is applied 
long after the microwave signal is detected. 



III. THEORY DETAILS 

We now give more details on the theory calculations. 
All these calculations are performed in the Holstein- 
Primakoff approximation, in which the spins and the 
resonator are described by harmonic oscillators, as ex- 
plained in the text. The system Hamiltonian is H/h = 
uj^{^)a^a + J2^ibjbj+J2'i9i{b\a — bja^), being the cou- 
pling constant of spin j with the resonator. 



A. Rabi oscillations 

We flrst calculate the oscillating signal resulting from 
the cycles of storage and retrieval of the resonator fleld 
of amplitude a into the spin ensemble (flg. 2d in the 
article). A flrst remark is that the system consisting of 
coupled harmonic oscillators, its time evolution does not 
depend on the initial fleld amplitude. As a result we 
only calculate aaabi = (0| a{t)a^{0) |0) which represents 
the probability amplitude that a photon created at t = 
is still present at time t. As shown in [2_ this quantity can 
be calculated by considering an effective non-Hermitian 
Hamiltonian 
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with complex angular frequencies Coq — uJr — and 
LUk — (^k ~ *7o/2 ; here, 70 is the spontaneous emission 
rate of each spin (that we take here to be zero since NV 
centers at low temperature have negligible energy relax- 
ation). Indeed, introducing the vector X{t) of coordi- 
nates [(a(t)at(0)),...,(6j(t)a^(0)),...] it can be shown 
that dX/dt — —{i/h)HeffX. The formal solution to this 
equation is then 



X{t)^C-'[is + ^H,ff/h)-'xm, 
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Figure 5: Measurement setup and wiring, (a) Full configuration of the measurement apparatus at room temperature, (b) 
Schematic of the wiring inside the dilusion refrigerator. LPFl, LPF2, LPF3, and LPF4 are low-pass filters having cutoff 
frequencies 1.35 GHz, 450 MHz, 450 MHz, and 1.35 GHz, respectively. BPFl actually consists of two series of low-pass (5.4 
GHz) and high-pass (1.3 GHz) filters, and BPF2 is a band-pass filter (2.5 to 4 GHz). CuNi coaxl, 2, and 4 are coaxial cables 
made of CuNi, and CuNi coax3 is a silver-plated CuNi coaxial cable. SC coax is a superconducting NbTi coaxial cable. Flex 
coax is a low-loss fiexible coaxial cable. The two circulators are STE-1502KS (Pamtech). Here the small rectangles represent 
ports terminated by 50 Q. HEMT is a cryogenic microwave amplifier having low-noise temprature (Caltech). A DC current for 
the DC magnetic field Bnv is applied through an RC low-pass filter at room temperature. 



with X{0) = xg and xq ^ (1,0, ...,0) . This 
implies that a{t) — xq^ ■ X{t) — Cr^ with 
ti(s) = xg"^ ■ (s -I- iHejfY^ ■ XG and C[f{s)] = 
/ f{t)dt, s being a complex number. Since ti(s) 
is not singular on its imaginary axis, we only need 
ti for pure imaginary argument s = —iuj to per- 
fom the Laplace transform inversion. As shown in 
[2], we have ii(— iw) = i/ [uj — ujq + injl — W{uj)] with 
W{uj) = gl^g j' p{uj')dLo' / [uj - bj' + i-fQ/2]. Computing 
a{t) is thus achieved by evaluating ti for the distribution 
p{lij) described in the main text, and numerically eva- 
luting the inverse Laplace transform. This calculation 
also enables us to calculate the transmission as plotted 
in Fig.le, because as shown in [2], the resonator trans- 
mission is obtained as tiuj) = — (k/2) ti{—iuj). 

An additional complication in the experiment is that 
the resonator-spin detuning 5 during the microwave pulse 
is finite, implying that the initial state is actually al- 
ready hybridized. To take this into account in the 
calculation, we make the approximation that this ini- 
tial state is a coherent superposition of Xg and xs = 
(0, gi, ■■,gN)/gens- The vector Xs is associated with 



the excitation of the superradiant mode. Thus the initial 
state is written x(t = 0) = cos{9/2)xg + ism{6 /2)xs 
with mixing angle tan6 = 2gens/S. Using a similar 
analysis, we now have anabi — cos(0/2)£~^ [ti{s)] + 
ism{e/2)C-'^ [t4{s)]. In addition to ti{s) we then need 
^4(5) = Xq{s + iHef f)~^xs, as shown in t4{~iuj) — 

-iti{-iuj)W{LL>) / gens- 



B. Ramsey fringes 

Each 7r/2 pulse necessary to realize the Ramsey-like ex- 
periment is realised by bringing the resonator and spins 
to resonance. For a fast pulse, the resonant interaction 
maps continuously the coherent state of the field to the 
superradiant mode in the spins. We calibrate the inter- 
action time in such a way as to transform the state xg 
into the superposition After the first 7r/2 pulse, 

the resonator is kept detuned from the spin ensemble for 
a time t. The system state at this point can be evalu- 
ated using eq. (2). We define Xcit) (resp. Xs{t)) as the 
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Figure 6: (a) Photograph of the area where the 4-SQUIDs 
array is placed in the resonator. Note that this photograph 
was taken before the SQUIDs were fabricated. The arrows 
shows a fiow of an applied fast current pulse, (b) Schematic 
of inside the dashed box in (a) (here the SQUIDs are drawn). 
The applied (return) current pulses is represented by a solid 
(dashed) arrow(s) . (c) Actual pulse sequence for the measure- 
ment: (i) excitation, (ii) interaction, (iii) data accumulation, 
and (iv) compensation in the sequence. After data sampling 
(iii), an inverted pulse (iv) is applied to compensate the time 
integral to (see text). The whole length of one sequence is 
5 us. 



C. Density distribution 

The spins density distribution that we use is the con- 
volution of a normalized Gaussian of standard deviation 
a with a normalized Lorentzian of HWHM 7. The re- 
sulting distribution is known as the Voigt profile and is 
also normalized. 



p{u}'; cr, 7) = / G{u}"; a)L{uj' - uj"; j)du}' 

.'2..,. (4) 



1 7 



00 Uj"'^/2<7'^ 



duj' . 



An important property is that we can compute analyt- 
ically the function W{uj) : 



W^(w;cr, 7, 70) = 



ay^TTj J {uj' — uj"y + ^'^ U} — uj' + i^Q 



(5) 



Taking the limit 70 — >■ we have 



vector of coordinates [(a(t)a^(O)) , (6j(i)a^(0)) , ...] at 
time t with initial conditions xq (resp. xs)- A second 
7r/2 pulse is then applied before the amplitude anp of 
the field in the resonator is measured : 

aRF - ^4 • U^niXait) - Xs{t)) 

^lixl + xl)-iXGit)-Xsit)) (3) 

= ^C-\h{s) - t2{s) + his) - t4{s)) , 

where ^2(5) — ■ (s -I- iH^ff)^^ ■ xs and t2{s) — xs^ ■ 
{s + iHeff)~^ -xq, ti and are defined above. As shown 
in [2, t2{-iuj) = -ti{-iuj)W{uj){s + iujoj/gl^^ and ts = 



^^"'"'^)=-^/2 

(6) 

The function that is taken in the theory plots of the 
main text actually takes into account the well-known NV 
center hyperfine splitting due to ^^iV by adding up three 
identical distributions p{ijj) separated by 2.2 MHz [3]. 
The resulting distribution is the one shown in figure 2 of 
the main text. This is conviently done since W[pi + P2] — 
W[pi\ + W[p2\ for any two different spin distributions 
Pi 2- Using this formula for W{lxj) we can evaluate the 
functions ti necessary for the simulation of the Rabi and 
Ramsey- like experiments. 



[1] Here a DC constant voltage is applied to the Q-port of the [2] I. Diniz et al, arxiv:quant-ph/1101.1842 (2011). 

IQ-mixer for the perpose of cancelling the offset, i.e., to [3] V. Acosta et al, Phys. Rev. B 80, 115202 (2009). 
prevent the microwave signal from leaking out during the 
pulse is off. 



